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0\ • Abstract 

(N 

C^ I We consider the full Dicke spin-boson model composed by a single bosonic mode and an en- 

semble of A^ identical two-level atoms with different couplings for the resonant and anti-resonant 
interaction terms, and incorporate a dipole-dipole interaction between the atoms. Assuming that 
the system is in thermal equilibrium with a reservoir at temperature /3~^, we compute the free 
. energy in the thermodynamic limit A^ — )■ oo in the saddle-point approximation to the path in- 
rN I tegral and determine the critical temperature for the superradiant phase transition. In the zero 
c^ • temperature limit, we recover the critical coupling of the quantum phase transition, presented in 
the literature. 
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1 Introduction 



At low temperatures a gas of N slow moving molecules with a few accessible electronic states 
interacting with the vacuum electromagnetic field in a cavity may be considered as a single quan- 
tum mechanical system with a ground state which determines the macroscopic properties of the 
gas and a set of excited states emerging as collective matter-field modes. In conditions where 
the dependence on the center of mass coordinates of the molecules may be omitted and for the 
simplest case of two level atoms, the system may be approximately described by Dicke [1] and 
generalized Dicke models which are written in terms of equivalent spin variables for the atoms. 
At zero temperature the properties of these models can be investigated by solving the quantum 
mechanical problem in a suitable approximation. In the original Dicke model [1] it is furthermore 
assumed that the atoms are confined to a small volume and that their direct interaction may be 
disregarded. The field is coupled to the total equivalent spin operator and enhanced rates for 
spontaneous radiation may be computed in the dipolar approximation for particular, so called 
superradiant states. In a modern view this is recognized as a hint of the existence of a quantum 
phase transition for this system [2] . A direct determination of this transition was obtained in Ref . 
[3] performing the diagonalization of the bimodal effective Hamiltonian resulting of applying a 
Holstein-Primakoff transformation [4, 5] in the thermodynamic limit and computing the critical 
coupling for which the boson expectation value vanish. In this work the energy of the excited 
collective modes was also computed and the relation of the quantum phase transition and the 
entanglement properties of the collective states was discussed. This connection is particularly 
interesting due to the wide interest that have attracted the applications to communications and 
quantum computing of entangled systems [6]. In this context efforts are focused in how to create 
and also how to detect entangled states in many-body systems [7, 8, 9], but many authors have also 
investigated the relation between entanglement of collective modes and quantum phase transitions 
[10, 11, 12, 13]. 

An alternative approach to the study of this model, which in fact was pursued much earlier is 
to explore the thermodynamical properties of the system at finite temperature by computing the 
partition function and the free energy. This was done first in the rotating wave approximation by 
Hepp and Lieb [14, 15] and Wang and Hioe [16] who determined the critical temperature for the 
superradiant phase transition in the Dicke model. In Ref. [17] a generalized Dicke model which 
include both the resonant and anti-resonant interaction terms was considered . The results in this 
works were confirmed in [18, 19] using different methods. 

The use of the path integral approach for the computation of the thermodynamical properties 
of these models was advocated by Popov and their collaborators [20, 21, 22, 23, 24] and produce 
results which are consistent with those obtained by other techniques in a very efficient way. The 
path integral formalism provides in particular a simple method to compute the spectrum of the 
system under consideration in either of it phases using the semiclassical expansion of the partition 
function [22, 25]. In Refs. [26, 27] these methods where used to determine the partition function 
and critical temperature of the generalized Dicke model and of the full Dicke model which emerges 



when different couplings for the resonant and anti-resonant terms are introduced. Then in Ref. 
[28] the spectrum of these models was determined by computing the quadratic approximation of 
the path integral. The zero temperature limit of the energies, in the case when the resonant and 
anti-resonant couplings are the same, were shown to correspond to the result in [3]. 

When dipole-dipole interactions are included the procedure using the Holstein Primakoff trans- 
formation is still useful and the critical coupling corresponding to the zero temperature quantum 
phase transition as well as the energy of the excited states may be computed [29, 30]. See also 
Refs. [31, 32, 33]. 

In this paper we compute the free energy in this generalized situation using functional methods 
and determine the critical temperature. Our result show a dependence of the critical coupling on 
the dipole-dipole coupling which in the zero temperature limit reproduce consistently the quantum 
phase transition point computed in [29]. The organization of the paper is as follows: In section II 
we discuss the interacting Hamiltonian of the model. In section III the functional integral for the 
model is performed and the critical temperature is determined. In section IV we present the free 
energy of the model. Finally, section V contains our conclusions. We use ks = c = h = 1 



2 The Hamiltonian and the fermionic model 

We consider the full Dicke model with dipole-dipole interaction. This model contemplates resonant 
and anti-resonant coupling terms between the atoms and the bosonic field with different couplings 
gi and g2. The Hamiltonian for this system, with the dipole-dipole interaction included reads 

X AT Q N 

ijtj 1=1 



In the above equation the operators a], af and erf = a^ satisfy the commutation relations [af , af] = 
2 e^'^'^ al for p,q,r = 1, 2, 3 and we use through the paper the operators af = | {aj ± i af). These 

operators satisfy [(Tj^,(T~] = af^^ and [(Tf,(Tf] = ±2o"f. The boson annihilation and creation 

operators h and 6''' satisfy the usual commutation relation rules. Each two-level atom interact with 
all other atoms of the ensemble with the same coupling strength \/N , and the summation is over 
all the atoms. Note that this model is related to models describing spin systems with long range 
interaction [34, 35, 36]. 

To proceed, let us define an auxiliary model to be called the fermionic full Dicke model in terms 
of fermionic raising and lowering operators aj, Ui^ (3j and /3j, that satisfy the anti-commutator 
relations ttjO;] + ajaj = 6ij and /3j/3J + /3J/3j = 6ij. We can also define the following bilinear 



combination of fermionic operators, alai — (3l(3i, al(5i and /3jai which obey the same commutation 
relations as the pseudo-spin operators (T?^-), at^^ and ct^n with the correspondence given by, 

af — > a\ai - /3J/3i , (2) 



a. 



+ . ^,t 



O",- 



ajlSi , (3) 

f3la,. (4) 

Substituting these relations in the Hamiltonian (1) we obtain the Hamiltonian Hp of the auxiliary 
fermionic full Dicke model with dipole-dipole coupling, 

\ N Q N 

ij i=l 

Af ^ 



The Hamiltonians if and Hp are defined in different Hilbert spaces, since each operator af in H 
acts on a two-dimensional Hilbert space, and each of the fermionic operators a], ai, (3j and /3j, 
in Hp acts on a four- dimensional Fock space. Nevertheless the partition function Z of the Dicke 
model may be written as a trace in the space of the fermionic operators due to the the following 
useful relation whose demonstration follows the one given by Popov for the Dicke model [21] 



Z = Tr[exp{-l3H)) =i'^Triexpi-l3HF iV^jj . (6) 

In this relation H is given by Eq. (1), Hp is given by Eq. (5) and the operator Np is defined by 

N 

Np = J2iala, + PJP,). (7) 

j=i 

It should be clear that the traces in Eq. (6) for each Hamiltonian are taken over their respective 
spaces. 



3 The partition function in the path integral approach 

As mentioned in the introduction the critical temperature for the Dicke and the full Dicke models 
may be computed using functional methods [27, 28]. In this section we consider the corresponding 
computation for the full Dicke model with dipole-dipole interaction. This problem was addressed 



earlier in Ref. [37], but within the approximation scheme proposed there no variation of the 
critical temperature was detected. In this section we prove that the introduction of the dipole- 
dipole interaction indeed modify the temperature for transition from normal to superradiant phase. 
Concerning the path integral approach that we use, we note that in this case due to the dipole- 
dipole interaction term the path integral is non Gaussian. To overcome this difficulty we introduce 
an auxiliary field in the path integral and perform a series expansion.. 

At finite temperature the Euclidean action S associated to the fermionic Dicke model is given 
by 

S = f^ dr U\t) drb{r) + E(«K^) dra,{T) + /3;(r) dMr))\ - j^ drHpir) , (8) 

where the Hamiltonian density Hf{t) is obtained from Eq. (5). In this case it takes the form 

Hf{t) = a;o6*(r)6(r) + -5:(^a*(r)a.(r) -/3;(r)/3,(r)j + 



N 

■'7= J 

N 



+ 



E«*(r)/3.(r)/3;(r)a,(r) + ^ 5: (a*(r) /3,(r) 6(r) + /3;(r) a,(r) 6*(r) ) + 
= E (/3;(r)a.(r)6(r) + a*(r) /3,(r) 6*(r)j. (9) 



Now we consider the formal ratio of the partition function of the dipole-dipole full Dicke model 
and the partition function of the free Dicke model. 



Z _ I[dr]ia,P,b)]exp[S-§!/^n{T)dr 



Zo J[di]{a, /3, b)] exp {So-§ j^ n{r)dr 
where the function n{T) in the above expression is defined by 



(10) 



N 

n{r) = E(<(^) ««(^) + A*(r)A(r)) . (11) 



i=l 



In Eq. (10) S = S'(6, 6*, a, q;^,/3, /3^) is the Euclidean action of the full Dicke model given by Eq. 
(8), 5*0 = 5*0 (&, b*, a, a^, (3, (3'^) is the free Euclidean action obtained by taking X = Qi = g2 = and 
[dri{a, I3,b)] is the functional measure. The functional integrals in Eq. (10), have to be done in 
the space of complex functions &*(r) and 6(r) and Grassmann variables a*{T), Q;j(r), /3*(t) and 
/3i{T). In thermal equilibrium in the imaginary time formalism, the integration variables in Eq. 
(10) obey periodic boundary conditions for the Bose field, i. e., 6(/3) = 6(0) and anti-periodic 
boundary conditions for Grassmann variables, i. e., q;j(/3) = — q;j(0) and /3j(/3) = — /3j(0). 



To advance in the computation of the quotient (10), we perform now the following decoupling 
transformation 

a,(r)^e^"a,(r), a*(r) ^ e" ^" a*(r) , 
A(r)-^e^^A(r), /3;(r) ^ e" ^^ /3*(r) . 

which makes the coefficient of the term proportional to ^(t) in the action to vanish. We should 
remark that the price to pay for this benefit is that the boundary conditions are modified. In Eq. 
(14), the Bose field still obeys periodic boundary conditions, i. e., b{/3) = 6(0), but the Fermi fields 
obey now the following boundary conditions: 

A(/3) = ^A(0), P:{P) = -tl3:{0). 

After this procedure we obtain 

Z_ ^ J[dv{a,(3,b)]e' 

Z, /[rfr7(a,/3,6)]e^o- ^ > 

In order to obtain the effective action of the bosonic mode we must integrate over the Grassmann 
Fermi fields in Eq. (14). The problem that we now are confronting is that in the action given by 
Eq. (9) there is a quartic term corresponding to the dipole-dipole interaction and the integral is 
non Gaussian. Although it is not possible to integrate this term directly, we can use an auxiliary 
field to circumvent this difficult and express the path integral as a series expansion. Since terms 
quadratic in the Grassman variables obey Bose statistics we have (Q;i(r) /3*(r)) (aj(r) /3*(r)) = 
{ajir) /3*{t)) (q;j(t) /3*(r)), and we can use following functional identity 

= No /'[ci^(r)]e-^°''K'*^"^'^"^ + ^^-^'^"^^'*^"^"'^"^ + ^^-i'*^"^"^*^"^^'^"0, (15) 

where [dr]{r)] is the functional measure for the functional integral with respect the fields r{T) and 

r*{T). The normalization factor is defined by Nq^ = J[dr]{r)] exp f /q dTr*{T) r(r) j, and the fields 

rir) and r*{T) satisfy the boundary conditions r(0) = r(/3) and r*(0) = r*(/3). Substituting this 
last expression in Z, we obtain 

Z = iVo 1^7^(0] [rfr7(a,/3, 6)] e^". (16) 

The new action Sr can be separated into a free action So{b,r) for the bosons and the field r{T) 
and a Gaussian fermionic part in the form, 

Sr = So{b,r)+ dTY.pl{T)Mr{b*,b)p,{T), (17) 



with 

So{b,r) = jyr {b*{r){dr-uJo)b{r) + r*{r)r{r)) , 

In (17), the column matrix Pi{T) arranges the Grassmann Fermi fields in the form 



P\ir) = (/3:(r) a*{r)) (19) 

and the matrix Mr{b*, b) is given in terms of the operators L = dr + fi/2 and L^, = dr — fi/2 by 



^ L -N-^/^(gib*{T)+g2b{T)-y/Xr{T)]^ 

^ -iV-i/2 (^g^ b (r) + g, b* (r) - VA r*(r) 



(20) 



From the form of Sr we see that the fermionic part of the functional integral given by Eq. (16) is 
Gaussian and can be done directly. Therefore we obtain 



N 



Z = No [dr]{r)] [dr]{b)] e^"(^'") (det M,(6*, b)) , (21) 

where in this case, [dri{b)] is the functional measure only for the bosonic field. With the help of 
the following property for matrices with operator components 

det I I = det [AD - ACA-^b) , (22) 

and using also determinant properties, we finally obtain, 

det M, (6*, 6) = 

det(LL,) det (l - N-^L;U gib + g2b* - \^r*)L-Ugib* + g^b - y/Xr)) . (23) 



Substituting Eq. (21) and Eq. (23) in Eq. (14), we have 

A = ^ (2A) 

Zo J[dv{r)] [dr]{b)] e^o(fe,r) ' ^ ' 



with Za defined by 

Za = J[dv{r)] [dr^m 

exp (So{b, r) + Ntr\n(l- N'^L;^ U b + g2b* - VXr*)L-^ (gi h* + g^h - TXr ) ) ) .(25) 



We have arrived at an effective action in terms of two bosonic modes interacting in a highly non 
trivial way. Since we are interested in knowing the asymptotic behavior of the quotient -^ in 
the thermodynamic limit, i. e., iV — )■ oo we analyze the asymptotic behavior of the last defined 
expression Za- First, to make the dependence in A^ more explicit, we rescale the bosonic fields by 

h -^ y/Nb, b* -^ VNb*, r -> \^ r and r* -^ ^/Nr*. We get 

Za = A{N) hd'qir)] [d'qih)] exp [n $(6, 6*, r, r*)) , (26) 

with the function $(6, 6*,r, r*) defined by 

$(6, 6*, r, r*) = So{b, r) + tr Inf 1 - L'Ugi b + g2b* - \^r*)L-^ (gi b* + g^b - ^/\r)\ . (27) 



The term A{N) in Eq. (26) comes from transforming the functional measures [dri{r)\ and [dri{b)\ 
under scaling the bosonic field by 6 — )> yNb^ b* — )■ ^/Nb*, r — )■ yN r and r* — )■ yN r* and will 
be reabsorbed when we transform back to the original variable. 

The asymptotic behavior of the functional integral appearing in Eq. (26) when A^ — )> oo, is 
obtained using the method of steepest descent [38]. We expand the function $(6*,6, r*,r) around 
stationary points 6(r) = b^ij), b*{T) 
points are solution of the equations 



6Q(r), rir) = ro(r) and r*{T) = rQ^r). The stationary 



S^{b,b*,r,r*) 



5b{T) 

5^{b,b*,r,r* 



bo,b^,ro,ro 







5^{b,b*,r,r* 



6r{T) 



6b*{T) 
5^{b,b*,r,r*] 



bo,b^,ro,r^ 







6r*{T) 



{2t 



which for the system at hand are constant functions, 6(r) = b^, &*(r) = b^, r{T) = Tq and r*{T) 
that fulfill the following system of algebraic equations, 

^o&o = (9i{9ibo + g2bl- y/Xr-Q^ + g2{gib*o + g2bo - ^/Xro)] — tanh l^^Aj , 



(29) 



'^o bl = (gi[gi bl + 5-2 &o - v^ro) + g2(gi &o + 5-2 &o - v^^o) ) 7^ tanh 



'P 



fi/ 



(30) 



v^ {91 K + g2bo- v^ro) — tanh | -^a j , 



ro = V\ (gi Bq + 5-2 &o - v^^o) 7;- tanh -^a , (31) 



r, 







VA [gi bo + g2K- yX r*) ^ tanh (^nA , (32) 

with Ha = \/^^ + 4: \gi bo + ^2 ^S - v^^o P- Substituting Eqs. (31) and (32) in Eqs. (29) and 
(30) we have also 

y/Xuobo = gir* + g2ro, 

V\uJob*Q = giro + g2r*o- (33) 

Substituting Eq. (33) in the Eq. (29) and imposing 60 7^ we finally obtain that 

^tanhf^^AV (34) 



(fi'i + fi'2)^ - ^0 A VLa V2 

This last equation allows us to calculate |6o|- Substituting Eq. (34) in Eq. (29) shows that 60 
is real and in consequence Tq is real. The critical temperature is characterized by the condition 
60 = 0. Therefore we have 

^"^ tanhf^f^V (35) 



(fl-i + g2Y - Wo A \2 

At zero temperature this result reproduces the critical coupling of the quantum phase transition 
discussed in Ref. [29]. We stress that the introduction of the dipole-dipole interaction modifies 
the critical temperature of the transition from the fluorescent to the superradiant phase. 



4 The free energy 

Let us complete the computation of the asymptotic behavior of the functional integral, appearing 
in Eq. (26), in the thermodynamic limit, A^ — )■ 00. We consider the two first leading terms in the 
functional integral appearing in Eq. (26) coming from the expansion of $(6, 6*,r, r*) around the 



maximal value 60, ^oj ''"0 and Tq this expansion is given by 

1 fl^ 
$(6, b*,r, r*) = $(60, bl, Tq, r^) + - dn dr^ 

1 Jo 

(bin) - bo , b*{n) - bl , r{n) - r^ , r*{n) - r*) 



Mi 



I b{r,)-bo \ 
h*{r2)-bl 
r{T2) - To 

V r*{T2)-r* J 



the matrix M$, is given by 

Sb{Ti)Sb(T2) 

5'^<S>{b,b* ,r,r*) 

T\/T _ Sb* (ti) Sb{T2) 

<5r(ri) Sb{T2) 
5^<J>(b,6*,r,r*) 



5^'!>{b,b*,r,r*) 

Sb{Tl) Sb* iT2) 

5'^^(b,b* ,r,r*) 

5b*(Tl)5b*(T2) 

5'^^(b,b* ,r,r*) 
<5r(Ti) 5b* (t2) 
5^<S>(b,b*,r,r*) 



5^'!>{b,b*,r,r*) 
56(ri) 5r{T2) 
P<S>{b,b*,r,r*) 
5b* (n) 5r(T2) 
5'^^(b,b* ,r,r*) 
5r(ri) 5r(r2) 
5'^^{b,b* ,r,r*) 



5'^<P{b,b*,r,r*) \ 
5fe(ri) 5r*(T2) 
5'^^(b,b* ,r,r*) 
56* (ti) 5r*{T2) 
5^<J>(6,fe*,r,r*) 
5r(ri) 5r*(T2) 
5^$(&,fe*,r,r*) 



\ 5r*(ri) 56(r2) 5r* (ti) Sb* (t2) 5r*(Ti) 5r(r2) 5r*(Ti) 5r*(T2) / 

Substituting this expansion given by Eq. (36) in Eq. (26) we have that 
Za = e^'^'^'oM [[dr]{r)] [dr]{b)] 
I 

1 tfi 

2 Jo 



bQ,b^,ro,r^ 



exp 



dn dT2 (bin) , b*{n) , r{n) , r*in)) M, 



(36) 



(37) 



(3^ 



r(r2) 
V \r*{r2))) 

To obtain the last expression, we applied the transformation 6(r) — 7'f6(r) + boj/y/N, b*(T) — )■ 

(6*(r) + blJ/VN, r{T) ^{r{T) + ro)/VN and r*{T) ^(r*{T) + r^)/VN to the fields in the 
functional integral. Terms of higher order in the fields in the expansion given by Eq. (36) generates 
terms in Eq. (38) suppressed by powers of 1/^/N. Therefore in the thermodynamic limit the Eq. 
(36) is a good approximation. 

The free energy F = —j^ In Z, in the thermodynamic limit, A^ — ?■ 00 is obtained from Eq. (24) 
using Eq. (38). In the normal phase F — Fo = 0, where Fq is the free energy for the non-interacting 
model. Finally, in the superradiant phase we have that 



F-Fn 



Wo (^A - ^^ 



4 ((^1 + ^2)' -A Wo) 



^■" 



'cosh(^^ 
cosh (^' 



(39) 



5 Conclusions 

Functional methods provide an efficient tool to compute the thermodynamics functions and critical 
properties of spin models. In the limit of zero temperature the results obtained by this means may 
be compared with the ones obtained directly from an exact or approximate quantum mechanical 
solution of the system involved. 

In this paper we compute the free energy of the full Dicke model with dipole-dipole interaction 
at temperature (3~^ and determine the critical temperature for the superradiant phase transition. 
To handle the non-Gaussian terms resulting from the presence of the dipole-dipole interaction we 
introduce a single auxiliary field which allow to compute an effective action in terms of two bosonic 
fields. This action showed to be manageable using the saddle point approximation and in this 
way we determine the corrections to the critical temperature of the full Dicke model computed 
in [17, 27]. This approach improves the approximation scheme of Ref. [37] where no correction 
was detected. For the full Dicke model with dipole-dipole interaction addressed in this paper 
the critical coupling associated with the quantum superradiant phase transition were obtained in 
Refs. [29, 30] following [3] by performing a Holstein Primakoff and then the diagonalization of the 
effective Hamiltonian so obtained. At zero temperature our result reproduces the critical coupling 
of the quantum phase transition discussed in Ref. [29, 30]. 

In Refs. [3, 29, 30] the relation of the entanglement of the collective modes of the Dicke 
models and the quantum phase transition was discussed. Since for finite temperature the thermal 
interaction usually induces decoherence in quantum systems it could be interesting to link this 
analysis with the thermodynamic behavior of the system. To this end, the excitation spectrum of 
the system should be computed and understood. This subject, as well as the effect of introducing 
disorder [39] in the dipole-dipole coupling are being investigated by the authors. 
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